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Abstract 



Sufficiently fast chemical equilibration of (strange) antibaryons in an en- 
vironment of nucleons, pions and kaons during the course of a relativistic 
heavy ion collision can be understood by a 'clustering' of mesons to buildt 
up baryon-antibaryon pairs. This multi-mesonic (fusion- type) process has to 
exist in medium due to the principle of detailed balance. Novel numerical 
calculations for a dynamical setup are presented. They show that - at maxi- 
mum SPS energies - yields of each antihyperon specie are obtained which are 
consistent with chemical saturated populations of T ~ 150 — 160 MeV, in line 
with popular chemical freeze-out parameters extracted from thermal model 
analyses. 
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I. BRIEF OVERVIEW ON ANTIHYPERON PRODUCTION 



Strangeness enhancement has been predicted a long time ago as a potential probe to find 
clear evidence for the temporary existence of a quark gluon plasma (QGP) in relativistic 
heavy ion collisions. A strong experimental effort has been made since and is still made for 
measuring strange particle abundancies in experiments at Brookhaven and at CERN (for 
a short recent review see [1]). In particular because of high production thresholds in bi- 
nary hadronic reaction channels antihyperons had been advocated as the appropriate QGP 
candidates [2]. Indeed, a satisfactory picture of nearly chemically saturated populations of 
antihyperons has been experimentally demonstrated over the last years with the Pb+Pb 
experiments NA49 and WA97 at CERN-SPS. For this statement, of course, a quantitative, 
theoretical analysis by employing a thermal (or 'statistical') model has to be invoked by 
fitting the thermodynamical parameters to the set of individual (strange) hadronic abun- 
dancies [3]. 

On the other hand, already since the first measurements with the lighter ions have been 
undertaken in the early nineties, the theoretical description of the antibaryon production 
within hadronic transport schemes in comparison to these data faced some severe difficulties. 
Some phenomenological motivated attempts to explain a more abundant production of anti- 
hyperons within a hadronic transport description [4] had been proposed like the appearance 
of color ropes, the fusion of strings, the percolation of strings, or the formation of high-dense 
hadronic clusters. The underlying mechanisms, however, have to be considered exotic and 
to some extent ad hoc, their purpose lies mainly to create (much) more antibaryon in the 
very early intial stage of the reaction (compared to simple rescaled p+p collisions). To some 
extent indeed the philosophy behind these mechanisms was to model a precursor of initial 
QGP formation within a hadronic transport scheme. On the other hand, in most of the 
transport calculations a dramatic role of subsequent antibaryon annihilation is observed, 
which, in return, has to be more than counterbalanced by these more exotic, initially occur- 
ing mechanisms. Large annihilation rates result if the free-space cross section is employed. 
The situation seems even more paradox with respect to the fact that the chemical descrip- 
tion within thermal models works indeed amazingly well for the antihyperons and being, of 
course, completely nondependent on the (large) magnitude of the annihilation cross section. 
For all of this reasons the theoretical and dynamical understanding of the production of 
(strange) antibaryons has remained a delicate and challenging task [5]. 
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FIGURES 




FIG. 1. Schematic picture for the multi-mesonic fusion-like reaction 3n + 2K *-* S + N. 

As we will demonstrate, a correct incorporation of the baryonic annihilation channels 
had actually not been done consistently. We have conjectured recently that a sufficiently 
fast redistributions of strange and light quarks into (strange) baryon-antibaryon pairs should 
be achieved by multi-mesonic fusion-type reactions of the type 

riiir + ri2K <-> Y + p (1) 

occuring in a moderately dense hadronic system [6] (illustrated in Fig. 1). The beauty of this 
argument lies in the fact that (at least) these special kind of multi-hadronic reactions have 
to be present because of the fundamental principle of detailed balance. As the annihilation 
of antihyperons on baryons is of dramatic relevance, the multi-mesonic (fusion-like) 'back- 
reactions' involving n\ pions and n<i kaons, where n 2 counts the number of anti-strange quarks 
within the antihyperon Y), must, in principle, be taken care of in a dynamical simulation. 
These reactions are then the most dominant source of production. This crucial fact had been 
overseen in all of the aforementioned treatments. The underlying reasoning was first raised 
by Rapp and Shuryak who described the maintenance of nearly perfect chemical equilibrium 
of antiprotons together with pions and nucleons during the very late stage of the expanding 
fireball before the particles loose contact [7]. 

The crucial input, though plausible, is now to assume that the annihilation cross sections 
for any strange or nonstrange antibaryon on any baryon are approximately the same order as 
for Np at the same relative momenta, i.e. &BY-^nw+n Y K ~ ^Np-^mr, being in the range of 50- 
100 mb for characteristic and moderately low momenta occuring in an expanding hadronic 
fireball. The equilibration timescale (r-p)( _1 ) ~ ^-/{^by v byPb) is to a good approximation 
proportional to the inverse of the density of baryons and their resonances. Adopting an initial 
density of 1-2 times normal nuclear matter density po f° r the initial and thermalized hadronic 
fireball, the antihyperons do equilibrate on a timescale of 1-3 fm/c well within the expansion 
timescale of the late hadronic fireball. Hence, fast chemical equilibration of the antihyperon 
abundancies is guaranteed by detailed balance with respect to the strong annihilation, the 
final yields then being independent of the actual size of the (large) annihilation cross section, 
solving the aforementioned paradox. 
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To be quantitative, (some) novel results by solving rate calculations for a dynamical setup 
are presented in the following. Before turning to their discussion, we will review briefly on 
some general ideas of the baryon-antibaryon annihilation process, in order to strengthen the 
one main assumption concerning the general size of its cross section, and on how to come 
to general master equations. 

II. ANTIBARYON ANNIHILATION AND ITS EFFECTIVE DESCRIPTION BY 

A MASTER EQUATION 

Antinucleon-nucleon annihilation is the strongest of all strong interaction processes. The 
strong annihilation of a nucleon and an antinucleon can be thought quantum mechanically 
as a complete absorbtive scattering process given approximately by the black disk formula 
&abs — ?r(-R + A) 2 or by a more sophisticated boundary condition description [8]. For p + p- 
annihilation one finds for the 'black disk' radius R — 1.07 fm, for which then this decription 
reproduces very accurately the total inelastic cross section as a function of the beam mo- 
mentum [8] and also its steep increase and diverging behaviour at low momenta. Hence, 
the quantum mechanical interpretation of the annihilation (being exotherm) is a picture of 
complete absorption. One is now tempted to generalize this simple and intuitive picture 
for all baryon-antibaryon annihilation processes. The only 'free' parameter which then 
can change is the radius R. Again, as it basically reflects the radius of the proton in the 
case of pp-annihilation, one would expect that this canonical value of 1 fm is of general 
validity for all the annihilation processes. Thus it is plausible to assume that the cross 
section for annihilation between any baryon and antibaryon should be rather the same for 
the same relative momenta. Indeed, there does exist old data on the total A + p cross sec- 
tion [9] with reasonably large value, although the A momenta was exceeding 4 GeV/c in all 
data taken. (These data are not conclusive to really extrapolate dwn to lower momenta, 
though they give a right indication: The authors of [9] reported an extrapolated formula of 
a abs = 47(±10)p A 1 ^ 2 mbGeV~ 1//2 , which is not fully analogue to the black disk formula, but 
has already a rather correct shape with a sizeable magnitude.) 

A more microscopic and quite popular picture, which can decribe quantitatively the 
complicated final states of individual outgoing mesonic channels, is the concept of two meson 
doorway states [10] given by 

B + B -> Mi + M 2 -* ni7r + n 2 K . (2) 

The main assumption is that the annihilation occurs exclusively via two mesons with near- 
est threshold dominance. Mi and M 2 can be rather highlying resonances and/or hybrid 
states which come close with their individual threshold. The final decay of these two mesons 
into the various channels of multiple pions and kaons is then treated statistically and mi- 
crocanonically. It turns out that such a description describes rather nicely reproduces the 
available data and can also be generalized straightforwardly for antihyperon annihilation 
[10]. There do exist a lot of other phenomenological and microscopic descriptions like quark 
models or flux tube models, which try to describe quantitatively the complicated process of 
the annihilation in more dynamical terms (for a review see [11]). 
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How can one incorporate such complicated processes into a description of transport 
dynamics? It is clear that as the very microscopy of the annihilation processes is not fully 
understood, one has to abandon any detailed local modelling, but has to turn to an effective 
coarse grained description which is guided by physical principles. Following the concepts of 
relativistic kinetic theory, the microscopic starting point is a Boltzmann-type equation of 
the form 

s , ff + i v/ > _ £ ^ J n j ^_ mt +PB E > 

{ni};B {rei,n 2 } {711,712} 

\(( ni ,n 2 \T\YB))\ 2 l(-)fyf B J] + I] - 

^ {"1,12} {«i,«2} J 

where 

{m} {111,12} \ {"1,12} / 

corresponds to the total annihilation cross section into the various possible multiple pion 
states. This form of the collision exactly incorporates the principle of detailed balance. 
The static fixed point of this equation (together with the usual binary kinetic processes) 
are thermal and chemical saturated distributions. Without the back reaction channels the 
equations would only have vanishing distributons as stable fixed points, which is, of course, 
unacceptable. Hence, the back reactions have to be considered as very basic and important 
ingredient of the transport description when trying to implement the baryon-antibaryon 
annihilation processes. 

Furthermore it is obvious that the back reactions will guarantee that the (strange) an- 
tibaryons become chemical saturated with the pions, kaons and nucleons on a very short 
timescale. To see this more explicit, and also for the further numerical treatment, one can 
bring the above Boltzmann equation into a more intuitive form of a master or rate equation. 
Assuming v re \ (Tyb(V^) to be roughly constant, which is actually a good approximation for 
the pp-annihilation, or, invoking a standard description of further effective coarse grain- 
ing by using thermally averaged cross sections and distributions, and furthermore taking 
the distributions in the Boltzmann approximantion, the following master equation for the 
respectively considered antihyperon density is obtained 



d Py = -{{(JybVyb)) \ PyPb - E M (ni<n2) (T, fi B , fi s )( P7T ) ni (p K ) n2 \ , (3) 



dt 



where the 'back-reactions' of several effectively clustering pions and kaons are incorporated 
in the 'mass-law' factor 

eg. eq. 

M (T \ — Y B 

m {n 1 ,n 2 )K- L ,P>B,Ps) — (p^-)™i(p^-)™2 Pm ' 

Here p ni (which will generally depend on the thermodynamical parameters) states the rela- 
tive probability of the reaction (1) to decay into a specific number ri\ of pions and ps denotes 
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the total number density of baryonic particles. As is well-known, the mass-law factor M 
depends only on the temperature and the baryon and strange quark chemical potentials. 
Ty = {(&yn v yn))Pb gives the effective annihilation rate of the respective antihyperon specie 
on a baryon. 
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FIG. 2. Time evolution of the (average) net baryon density for midrapidity |AY| < 1 and 
central Pb+Pb-collision at 160 AGeV obtained within a dynamical transport simulation. Here the 
amount of baryon number residing still in string-like excitations is explicitely discarded. String-like 
excitations have disappeared after about 3-4 fm/c, so that from this time on a pure hadronic fireball 
develops and expands. Its initial net baryon-density starts slightly above ps = %Po- 

For a further manipulation one has to make assumptions for the various abundancies 
occuring in the master equation. Nonequilibrium inelastic hadronic reactions can explain to 
a good extent the overall strangeness production seen experimentally: The major amount of 
the produced kaons at SPS-energies can be understood in terms of still early and energetic 
non-equilibrium interactions [12]. In Fig. 2 the time evolution of netbaryon density at midra- 
pidiy obtained within a dynamical transport simulation [13] is depicted. All strangeness is 
being produced still when string-like excitations are governing the dynamics. In the later 
hadronic stage the number of strange quarks stays more or less constant and can only be 
redistributed among the various hadrons [14]. Also the pions and nucleons do stay more 
or less at thermal equilibrium. Refering to the master equation (3), one can then take 
the pions, baryons and kaons to stay approximately in thermal equilibrium throughout the 
later hadronic evolution of the collision, the later being modelled to be an isentropic expan- 
sion with fixed total entropy content being specified via the entropy per baryon ratio S/A 
(compare with Fig. 3). (3) has the intuitive form 

j t Pv = -r Y { PY -Py 1 } . (4) 

The production rate per unit volume dN Y /dtdV is given by the non- vanishing, though small 
value r Y p^- Still this rate is enough to populate the antihyperons to their (very small) 
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equilibrium value. 
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FIG. 3. The netbaryon density pg of a chemically fully equilibrated thermal hadronic resonance 
gas as function of decreasing temperature for the situation of an isentropic expansion, i.e. for 
various constant entropy /baryon ratios. For the value S/A = 250 (a situation expected at RHIC) 
the total baryon number density p nuc is plotted instead. 

For the calculations solving this master equation (4) one has to employ an 'effective' 
volume V(t) in order to extrapolate from Fig. 2 to Fig. 3, i.e. simulating the global charac- 
teristic of the expansion and the dilution of the baryon density and thus the annihilation rate 
Ty(t). The 'effective' (global or at midrapidity) volume V(t) is parametrized as function of 
time by longitudinal Bjorken expansion and including a transversal expansion either with a 
linear profile 



(taking appropriate parameter to simulate slower or faster expansion) or with an 

accelerating radial flow 



Veff,lin(t > t ) 



n(ct) (R + v Un (t -t )) 2 



(5) 



V e ff,acc(t >t ) = 71 (ct) (Rq + v Q (t - t ) + 0.5a (t - t ) 2 ) 



(6) 



with R = 6.5 fm, v = 0.15 c and a = 0.05 c 2 



/ fm for the later modeling. 
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III. RESULTS AND IMPLICATIONS 
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FIG. 4. The anti-A to baryon number ratio N^/Nb (t) as a function of time for various velocity- 
parameters vu n for the transverse expansion. The entropy per baryon is taken as S/A = 30, to = 3 
fm/c and T = 190 MeV. 

At starting time t an initial temperature T is chosen. (T is set to 190 MeV for 
the SPS and 150MeV for the AGS situation, while the initial energy densities are then 
about 1 GeV/fm 3 .) From (5) or (6) together with the constraint of conserved entropy the 
temperature and the chemical potentials do follow as function of time and thus also pB{t) 
as well as p^{t) within the hadronic resonance gas. As a minimal assumption the initial 
abundancy of antihyperons is set to zero. Equation (4), taking into account the volume 
dilution, is solved for each specie. For the thermally averaged cross section we take a simple 
constant value of (a ann v) := a = 40 mb. 
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FIG. 5. The anti-A to baryon number ratio N^/Ng (t) as a function of time for various im- 
plemented annihiation cross section a e ff = Xctq. The entropy per baryon is taken as S/A = 30, 
to = 3 fm/c and T = 190 MeV. 

In Fig. 4 the number of As (normalized to the conserved net baryon number) as a function 
of time is depicted. The entropy per baryon is chosen as S/A = 30 being characteristic to 
global ('47r') SPS results [15]. Here we have chosen the ansatz (5)for the expansion of the 
volume being linear in time in the transverse direction. The parameter vu n is varied to 
simulate slow or fast expansion of the late hadronic fireball. The general characteristics 
is that first the antihyperons are dramatically being populated, and then in the very late 
expansion some more are still being annihilated, depending on how fast the expansion goes. 
A rapid expansion gives a higher yield, which can increase the final yield by a factor of 2 to 3. 
However, the typical expansion behaviour obtained from simulations or extracted from the 
analysis of transverse momentum slopes of individual hadrons (pions and protons) is that 
at the late stages the transverse expansion velocity shoul be about 0.5 c. In the following 
we stay to the second ansatz (6) which extrapolates from an initially slower to a later faster 
expansion. 
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FIG. 6. The antihyperon to baryon number ratio Ny/N B (T) and /N B (T) (dotted line) 
as a function of the decreasing temperature. Parameters are the same as in Fig. 5. 

In Fig. 5 the number of As as a function of time is depicted, where now the cross section 
employed is varied by a constant factor, i.e. a e ff = \cr . The results are rather robust 
against a variation by a factor of 2 in the cross section. Typically (for A = 1) about more 
than 5 times in number of antihyperons are created during the evolution compared to the 
final number freezing out, thus reflecting the fast back ('annihilation') and forth ('creation') 
processes at work dictated by detailed balance. 
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FIG. 7. N-^/Nb (t) as a function of time for various entropy content described via the entropy 
per baryon ratio (S/A = 20 — 40). Other parameters are as in Fig. 5. 

In Fig. 6 the number of antihyperons of each specie are now shown as a function of the 
decreasing temperature T(t) of the hadronic system. For a direct comparison the instanta- 
neous equilibrium abundancy N^'(T(t), /Ab(£), /i s (t))/NB is also given. As noted above, after 
a fast initial population, the individual yields of the antihyperons do overshoot their respec- 
tive equilibrium number and then do finally saturate at some slightly smaller value. More- 
over, one notices that the yields effectively do saturate at a number which can be compared 
to an equivalent equilibrium number at a temperature parameter around T e ff 150 — 160 
MeV, being strikingly close to the ones obtained within the various thermal analyses [3]. 

In Fig. 7 the number of anti-As as a function of time is given for various entropy per 
baryon ratios. One notices that the final value in the yield significantly depends on the 
entropy content, or, in other words, on the baryochemical potential. We note that the results 
at midrapidity from WA97 can best be reproduced by employing an entropy to baryon ratio 
S/A — 40. Indeed, at midrapidity one qualitatively expects a higher entropy content due 
to the larger pion to baryon ratio as compared to full £ 47r' data over all rapidities. At this 
point it will also be very interesting to compare our semi- quantitative calculations with the 
new results from NA49 on the A-yield at lower SPS energies of 80 AGeV and 40 AGeV 
with lower entropy contents, respectively. As the presented results are very sensitive on the 
entropy content, one first needs a clean analysis to obtain a rather accurate S/A number 
from the measured pion and proton abundancies. 
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FIG. 8. Nj^/Nb (T) and N^' /Nb(T) as a function of decreasing temperature for a characteris- 
tic AGS situation with an entropy content of S/A = 12 for various implemented annihiation cross 
section a e ff = \ao- *o = 5 fm/c and To = 150 MeV. 

There is also a clear hint at AGS energies of enhanced anti-A production [16]. For most 
central collisions more anti-As are found compared to anti-protons, which is quite puzzling 
as within a thermal model analysis this ratio is found not to be larger than 1. This enhanced 
ratio of anti-As compared to anti-protons at AGS energies one can understand in a way that 
one assumes that their annihilation cross section on baryons is just slightly smaller than for 
the antiprotons. In Fig. 8 a similar study like that of Fig. 5 is shown for a characteristic 
situation at AGS. For smaller, yet not too small effective cross sections the final yield can 
here be enlarged by a factor of 2 compared to the case with a 'full' crossection, as the final 
reabsorption is not as effective. But, of course, this idea is speculation at present. Also, 
we remark that the As effectively do saturate at an equivalent equilibrium number at a 
temperature parameter around T e ff xs 120 — 130 MeV. Unfortunately, there are no data for 
S at AGS. Again the new NA49 data at lower energies are worthwhile to pursue. A detailed 
measurement of all antihyperons represents also an excellent opportunity for future heavy 
ion facilities at an energy upgraded GSI. 

To summarize, multi-mesonic production of antihyperons is a consequence of detailed 
balance and, as the annihilation rate is large, it is by far the most dominant source in a 
hadronic gas. This is a remarkable observation, as it clearly demonstrates the importance 
of hadronic multi-particle channels, occuring frequently enough in a (moderately) dense 
hadronic environment in order to populate and chemically saturate the rare antibaryons. In 
order to be more competitive for a direct comparison with various experimental findings, new 
strategies have to be developed to describe for such multi-particle interactions within present 
day transport codes. A significant first step forward was very recently made; first results 
concerning the production of anti-protons at AGS and SPS energies are quite impressive 
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[17]. Another strategy could be to exploit microscopically the concept of two meson doorway 
states (2) and their sequential decay by standard binary scattering processes. 
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BRIEF OVERVIEW ON ANTIHYPERON PRODUCTION 

Strangeness enhancement has been predicted a long time ago as a potential probe to find 
clear evidence for the temporary existence of a quark gluon plasma (QGP) in relativistic 
heavy ion collisions. A strong experimental effort has been made since and is still made 
for measuring strange particle abundancies in experiments at Brookhaven and at CERN 
(for a short recent review see [1]). In particular because of high production thresholds in 
binary hadronic reaction channels antihyperons had been advocated as the appropriate 
QGP candidates [2]. Indeed, a satisfactory picture of nearly chemically saturated pop- 
ulations of antihyperons has been experimentally demonstrated over the last years with 
the Pb+Pb experiments NA49 and WA97 at CERN-SPS. For this statement, of course, a 
quantitative, theoretical analysis by employing a thermal (or 'statistical') model has to 
be invoked by fitting the thermodynamical parameters to the set of individual (strange) 
hadronic abundancies [3]. 

On the other hand, already since the first measurements with the lighter ions have been 
undertaken in the early nineties, the theoretical description of the antibaryon production 
within hadronic transport schemes in comparison to these data faced some severe dif- 
ficulties. Some phenomenological motivated attempts to explain a more abundant pro- 
duction of antihyperons within a hadronic transport description [4] had been proposed 
like the appearance of color ropes, the fusion of strings, the percolation of strings, or 
the formation of high-dense hadronic clusters. The underlying mechanisms, however, 
have to be considered exotic and to some extent ad hoc, their purpose lies mainly to 
create (much) more antibaryon in the very early intial stage of the reaction (compared 
to simple rescaled p+p collisions). To some extent indeed the philosophy behind these 
mechanisms was to model a precursor of initial QGP formation within a hadronic trans- 



FIGURE 1. Schematic picture for the multi-mesonic fusion-like reaction 3ji + 2K -h- S + N. 



port scheme. On the other hand, in most of the transport calculations a dramatic role of 
subsequent antibaryon annihilation is observed, which, in return, has to be more than 
counterbalanced by these more exotic, initially occuring mechanisms. Large annihila- 
tion rates result if the free-space cross section is employed. The situation seems even 
more paradox with respect to the fact that the chemical description within thermal mod- 
els works indeed amazingly well for the antihyperons and being, of course, completely 
nondependent on the (large) magnitude of the annihilation cross section. For all of this 
reasons the theoretical and dynamical understanding of the production of (strange) an- 
tibaryons has remained a delicate and challenging task [5]. 

As we will demonstrate, a correct incorporation of the baryonic annihilation channels 
had actually not been done consistently. We have conjectured recently that a sufficiently 
fast redistributions of strange and light quarks into (strange) baryon-antibaryon pairs 
should be achieved by multi-mesonic fusion-type reactions of the type 

n l n-\~n 2 K <H> Y + p (1) 

occuring in a moderately dense hadronic system [6] (illustrated in Fig. 1). The beauty of 
this argument lies in the fact that (at least) these special kind of multi-hadronic reactions 
have to be present because of the fundamental principle of detailed balance. As the 
annihilation of antihyperons on baryons is of dramatic relevance, the multi-mesonic 
(fusion-like) 'back-reactions' involving n l pions and n 2 kaons, where n 2 counts the 
number of anti- strange quarks within the antihyperon Y), must, in principle, be taken 
care of in a dynamical simulation. These reactions are then the most dominant source of 
production. This crucial fact had been overseen in all of the aforementioned treatments. 
The underlying reasoning was first raised by Rapp and Shuryak who described the 
maintenance of nearly perfect chemical equilibrium of antiprotons together with pions 
and nucleons during the very late stage of the expanding fireball before the particles 
loose contact [7]. 

The crucial input, though plausible, is now to assume that the annihilation cross 
sections for any strange or nonstrange antibaryon on any baryon are approximately the 
same order as for Np at the same relative momenta, i.e. a B f^. nn+n K ~ a Np^tac> being 
in the range of 50-100 mb for characteristic and moderately low momenta occuring in 
an expanding hadronic fireball. The equilibration timescale (r^)^ 1 ) ~ l/(o By v BY p B ) 
is to a good approximation proportional to the inverse of the density of baryons and 



their resonances. Adopting an initial density of 1-2 times normal nuclear matter density 
p for the initial and thermalized hadronic fireball, the antihyperons do equilibrate on a 
timescale of 1-3 fm/c well within the expansion timescale of the late hadronic fireball. 
Hence, fast chemical equilibration of the antihyperon abundancies is guaranteed by 
detailed balance with respect to the strong annihilation, the final yields then being 
independent of the actual size of the (large) annihilation cross section, solving the 
aforementioned paradox. 

To be quantitative, (some) novel results by solving rate calculations for a dynamical 
setup are presented in the following. Before turning to their discussion, we will review 
briefly on some general ideas of the baryon-antibaryon annihilation process, in order to 
strengthen the one main assumption concerning the general size of its cross section, and 
on how to come to general master equations. 



ANTIBARYON ANNIHILATION AND ITS EFFECTIVE 
DESCRIPTION BY A MASTER EQUATION 

Antinucleon-nucleon annihilation is the strongest of all strong interaction processes. The 
strong annihilation of a nucleon and an antinucleon can be thought quantum mechani- 
cally as a complete absorbtive scattering process given approximately by the black disk 
formula o abs = n(R + X) 2 or by a more sophisticated boundary condition description 
[8]. For p + p-annihilation one finds for the 'black disk' radius R = 1.07 fm, for which 
then this decription reproduces very accurately the total inelastic cross section as a func- 
tion of the beam momentum [8] and also its steep increase and diverging behaviour at 
low momenta. Hence, the quantum mechanical interpretation of the annihilation (being 
exotherm) is a picture of complete absorption. One is now tempted to generalize this 
simple and intuitive picture for all baryon-antibaryon annihilation processes. The only 
'free' parameter which then can change is the radius R. Again, as it basically reflects the 
radius of the proton in the case of ^-annihilation, one would expect that this canonical 
value of 1 fm is of general validity for all the annihilation processes. Thus it is plausi- 
ble to assume that the cross section for annihilation between any baryon and antibaryon 
should be rather the same for the same relative momenta. Indeed, there does exist old 
data on the total A + p cross section [9] with reasonably large value, although the A mo- 
menta was exceeding 4 GeV/c in all data taken. (These data are not conclusive to really 
extrapolate dwn to lower momenta, though they give a right indication: The authors of 
[9] reported an extrapolated formula of o abs = 47(±10)p~ 1,/2 mbGeV -1 / 2 , which is not 
fully analogue to the black disk formula, but has already a rather correct shape with a 
sizeable magnitude.) 

A more microscopic and quite popular picture, which can decribe quantitatively the 
complicated final states of individual outgoing mesonic channels, is the concept of two 
meson doorway states [10] given by 

B + B ->■ M l +M 2 ->■ n l K-\-n 2 K . (2) 

The main assumption is that the annihilation occurs exclusively via two mesons with 
nearest threshold dominance. M 1 and M 2 can be rather highlying resonances and/or 



hybrid states which come close with their individual threshold. The final decay of 
these two mesons into the various channels of multiple pions and kaons is then treated 
statistically and microcanonically. It turns out that such a description describes rather 
nicely reproduces the available data and can also be generalized straightforwardly for 
antihyperon annihilation [10]. There do exist a lot of other phenomenological and 
microscopic descriptions like quark models or flux tube models, which try to describe 
quantitatively the complicated process of the annihilation in more dynamical terms (for 
a review see [11]). 

How can one incorporate such complicated processes into a description of transport 
dynamics? It is clear that as the very microscopy of the annihilation processes is not 
fully understood, one has to abandon any detailed local modelling, but has to turn to an 
effective coarse grained description which is guided by physical principles. Following 
the concepts of relativistic kinetic theory, the microscopic starting point is a Boltzmann- 
type equation of the form 

Y {«[};£ F v ' O {n l ,n 2 } K ' i \ n \Pis 

\((n u n 2 \T\YB))\ 2 \(-)f ? f B ft (I + fi) + J] W ~ fyW ~ f B ) 



where 



tlUX {n l }{n v n 2 } J W lE i \ {n { ,n 2 } I 



YB 



corresponds to the total annihilation cross section into the various possible multiple pion 
states. This form of the collision exactly incorporates the principle of detailed balance. 
The static fixed point of this equation (together with the usual binary kinetic processes) 
are thermal and chemical saturated distributions. Without the back reaction channels the 
equations would only have vanishing distributons as stable fixed points, which is, of 
course, unacceptable. Hence, the back reactions have to be considered as very basic and 
important ingredient of the transport description when trying to implement the baryon- 
antibaryon annihilation processes. 

Furthermore it is obvious that the back reactions will guarantee that the (strange) 
antibaryons become chemical saturated with the pions, kaons and nucleons on a very 
short timescale. To see this more explicit, and also for the further numerical treatment, 
one can bring the above Boltzmann equation into a more intuitive form of a master or 
rate equation. Assuming v rd Oy B {\/s) to be roughly constant, which is actually a good 
approximation for the pp-annihilation, or, invoking a standard description of further 
effective coarse graining by using thermally averaged cross sections and distributions, 
and furthermore taking the distributions in the Boltzmann approximantion, the following 
master equation for the respectively considered antihyperon density is obtained 



J t P? = -^yb v yb)) { P?Pb ~ I M {ni n 2) (T,n B ,n s )( Pn ) n i (p K ) n 2 } , (3 ) 
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FIGURE 2. Time evolution of the (average) net baryon density for midrapidity | AY | < 1 and central 
Pb+Pb-collision at 160 AGeV obtained within a dynamical transport simulation. Here the amount of 
baryon number residing still in string-like excitations is explicitely discarded. String-like excitations have 
disappeared after about 3-4 fm/c, so that from this time on a pure hadronic fireball develops and expands. 
Its initial net baryon-density starts slightly above p B = 2p . 

where the 'back-reactions' of several effectively clustering pions and kaons are incorpo- 
rated in the 'mass-law' factor 

n eq. n eq. 

" \ Y B 

M {n v n 2 )y T ^ B iVs) = ^ p eq.y^ pe q.y 2 Pn, ■ 

Here p Ux (which will generally depend on the thermodynamical parameters) states the 
relative probability of the reaction (1) to decay into a specific number n 1 of pions and p B 
denotes the total number density of baryonic particles. As is well-known, the mass-law 
factor M depends only on the temperature and the baryon and strange quark chemical 
potentials. T ? = {{^yn v yn^Pb gi ves the effective annihilation rate of the respective 
antihyperon specie on a baryon. 

For a further manipulation one has to make assumptions for the various abundancies 
occuring in the master equation. Nonequilibrium inelastic hadronic reactions can explain 
to a good extent the overall strangeness production seen experimentally: The major 
amount of the produced kaons at SPS-energies can be understood in terms of still early 
and energetic non-equilibrium interactions [12]. In Fig. 2 the time evolution of netbaryon 
density at midrapidiy obtained within a dynamical transport simulation [13] is depicted. 
All strangeness is being produced still when string-like excitations are governing the 
dynamics. In the later hadronic stage the number of strange quarks stays more or less 
constant and can only be redistributed among the various hadrons [14]. Also the pions 
and nucleons do stay more or less at thermal equilibrium. Refering to the master equation 
(3), one can then take the pions, baryons and kaons to stay approximately in thermal 
equilibrium throughout the later hadronic evolution of the collision, the later being 
modelled to be an isentropic expansion with fixed total entropy content being specified 
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FIGURE 3. The netbaryon density p B of a chemically fully equilibrated thermal hadronic resonance gas 
as function of decreasing temperature for the situation of an isentropic expansion, i.e. for various constant 
entropy/baryon ratios. For the value S/A = 250 (a situation expected at RHIC) the total baryon number 
density p nuc is plotted instead. 

via the entropy per baryon ratio S/A (compare with Fig. 3). (3) has the intuitive form 

SPp = -Mp,-P?}. (4) 

The production rate per unit volume dN y /dtdV is given by the non- vanishing, though 
small value r ? pf 9 . Still this rate is enough to populate the antihyperons to their (very 
small) equilibrium value. 

For the calculations solving this master equation (4) one has to employ an 'effec- 
tive' volume V(t) in order to extrapolate from Fig. 2 to Fig. 3, i.e. simulating the global 
characteristic of the expansion and the dilution of the baryon density and thus the annihi- 
lation rate T Y (t). The 'effective' (global or at midrapidity) volume V(t) is parametrized 
as function of time by longitudinal Bjorken expansion and including a transversal ex- 
pansion either with a linear profile 

V effM {t > t Q ) = K{ct) (R + Vlin (t-t )) 2 (5) 

(taking v Un as an appropriate parameter to simulate slower or faster expansion) or with 
an accelerating radial flow 

^//,acc( ? >%) = X( ct ) (*0 + v o('-'o)+°-H('-'o) 2 ) 2 ( 6 ) 

with R = 6.5 fm, v Q = 0.15c and a Q = 0.05 c 2 / fm for the later modeling. 
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FIGURE 4. The anti-A to baryon number ratio N^/N B (t) as a function of time for various velocity 
parameters v [in for the transverse expansion. The entropy per baryon is taken as S/A = 30, t = 3 fm/c and 
/;, L90MeV. 

RESULTS AND IMPLICATIONS 

At starting time t Q an initial temperature T Q is chosen. (T is set to 190 MeV for the 
SPS and \50MeV for the AGS situation, while the initial energy densities are then 
about 1 GeV/fm 3 .) From (5) or (6) together with the constraint of conserved entropy 
the temperature and the chemical potentials do follow as function of time and thus also 
p B {t) as well as p e - q {t) within the hadronic resonance gas. As a minimal assumption 
the initial abundancy of antihyperons is set to zero. Equation (4), taking into account the 
volume dilution, is solved for each specie. For the thermally averaged cross section we 
take a simple constant value of ((7 ann v) := o Q = 40 mb. 

In Fig. 4 the number of As (normalized to the conserved net baryon number) as a 
function of time is depicted. The entropy per baryon is chosen as S/A = 30 being 
characteristic to global ('47r') SPS results [15]. Here we have chosen the ansatz (5)for the 
expansion of the volume being linear in time in the transverse direction. The parameter 
v Un is varied to simulate slow or fast expansion of the late hadronic fireball. The general 
characteristics is that first the antihyperons are dramatically being populated, and then 
in the very late expansion some more are still being annihilated, depending on how fast 
the expansion goes. A rapid expansion gives a higher yield, which can increase the final 
yield by a factor of 2 to 3. However, the typical expansion behaviour obtained from 
simulations or extracted from the analysis of transverse momentum slopes of individual 
hadrons (pions and protons) is that at the late stages the transverse expansion velocity 
shoul be about 0.5 c. In the following we stay to the second ansatz (6) which extrapolates 
from an initially slower to a later faster expansion. 

In Fig. 5 the number of As as a function of time is depicted, where now the cross 
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FIGURE 5. The anti-A to baryon number ratio N^/N B (t) as a function of time for various implemented 
annihiation cross section a e ^ = A a . The entropy per baryon is taken as S/A = 30, f = 3 fm/c and 
r = 190MeV. 
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FIGURE 6. The antihyperon to baryon number ratio N f /N B (T) and N e . q - /N B (T) (dotted line) as a 
function of the decreasing temperature. Parameters are the same as in Fig. 5. 
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FIGURE 7. N^/N B (t) as a function of time for various entropy content described via the entropy per 
baryon ratio (S/A = 20 — 40). Other parameters are as in Fig. 5. 

section employed is varied by a constant factor, i.e. o e ^ = X o Q . The results are rather 
robust against a variation by a factor of 2 in the cross section. Typically (for X = 1) about 
more than 5 times in number of antihyperons are created during the evolution compared 
to the final number freezing out, thus reflecting the fast back ('annihilation') and forth 
('creation') processes at work dictated by detailed balance. 

In Fig. 6 the number of antihyperons of each specie are now shown as a function 
of the decreasing temperature T(t) of the hadronic system. For a direct comparison 
the instantaneous equilibrium abundancy N^ q -(T(t),iJ, B (t) } iJ, s (t))/N B is also given. As 
noted above, after a fast initial population, the individual yields of the antihyperons 
do overshoot their respective equilibrium number and then do finally saturate at some 
slightly smaller value. Moreover, one notices that the yields effectively do saturate at a 
number which can be compared to an equivalent equilibrium number at a temperature 
parameter around T g ^ 150— 160 MeV, being strikingly close to the ones obtained 
within the various thermal analyses [3]. 

In Fig. 7 the number of anti-As as a function of time is given for various entropy per 
baryon ratios. One notices that the final value in the yield significantly depends on the 
entropy content, or, in other words, on the baryochemical potential. We note that the 
results at midrapidity from WA97 can best be reproduced by employing an entropy to 
baryon ratio S/A = 40. Indeed, at midrapidity one qualitatively expects a higher entropy 
content due to the larger pion to baryon ratio as compared to full '4tt' data over all 
rapidities. At this point it will also be very interesting to compare our semi-quantitative 
calculations with the new results from NA49 on the A-yield at lower SPS energies of 80 
AGeV and 40 AGeV with lower entropy contents, respectively. As the presented results 
are very sensitive on the entropy content, one first needs a clean analysis to obtain a 
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FIGURE 8. N^/N B (T) and N e ^ /N B {T) as a function of decreasing temperature for a characteristic 
AGS situation with an entropy content of S/A = 12 for various implemented annihiation cross section 
a eff = X a Q . t = 5 fm/c and T = 150 MeV. 

rather accurate S/A number from the measured pion and proton abundancies. 

There is also a clear hint at AGS energies of enhanced anti-A production [16]. For 
most central collisions more anti-As are found compared to anti-protons, which is quite 
puzzling as within a thermal model analysis this ratio is found not to be larger than 
1. This enhanced ratio of anti-As compared to anti-protons at AGS energies one can 
understand in a way that one assumes that their annihilation cross section on baryons 
is just slightly smaller than for the antiprotons. In Fig. 8 a similar study like that of 
Fig. 5 is shown for a characteristic situation at AGS. For smaller, yet not too small 
effective cross sections the final yield can here be enlarged by a factor of 2 compared 
to the case with a 'full' crossection, as the final reabsorption is not as effective. But, 
of course, this idea is speculation at present. Also, we remark that the As effectively 
do saturate at an equivalent equilibrium number at a temperature parameter around 
T £ rr m 120 — 130 MeV. Unfortunately, there are no data for S at AGS. Again the new 
NA49 data at lower energies are worthwhile to pursue. A detailed measurement of all 
antihyperons represents also an excellent opportunity for future heavy ion facilities at an 
energy upgraded GSI. 

To summarize, multi-mesonic production of antihyperons is a consequence of detailed 
balance and, as the annihilation rate is large, it is by far the most dominant source in a 
hadronic gas. This is a remarkable observation, as it clearly demonstrates the importance 
of hadronic multi-particle channels, occuring frequently enough in a (moderately) dense 
hadronic environment in order to populate and chemically saturate the rare antibaryons. 
In order to be more competitive for a direct comparison with various experimental find- 
ings, new strategies have to be developed to describe for such multi-particle interactions 



within present day transport codes. A significant first step forward was very recently 
made; first results concerning the production of anti-protons at AGS and SPS energies 
are quite impressive [17]. Another strategy could be to exploit microscopically the con- 
cept of two meson doorway states (2) and their sequential decay by standard binary 
scattering processes. 
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